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ABSTRACT 


Eukaryotic genome size data are important both as 
the basis for comparative research into genome 
evolution and as estimators of the cost and difficulty 
of genome sequencing programs for non-model 
organisms. In this study, the genome size of 14 
species of fireflies (Lampyridae) (two genera in 
Lampyrinae, three genera in Luciolinae, and one 
genus in subfamily incertae sedis) were estimated 
by propidium iodide (Pl)-based flow cytometry. The 
haploid genome sizes of Lampyridae ranged from 
0.42 to 1.31 pg, a 3.1-fold span. Genome sizes of 
the fireflies varied within the tested subfamilies and 
genera. Lamprigera and Pyrocoelia species had 
large and small genome sizes, respectively. No 
correlation was found between genome size and 
morphological traits such as body length, body width, 
eye width, and antennal length. Our data provide 
additional information on genome size estimation of 
the firefly family Lampyridae. Furthermore, this study 
will help clarify the cost and difficulty of genome 
sequencing programs for non-model organisms and 
will help promote studies on firefly genome evolution. 
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INTRODUCTION 


Fireflies, in the family Lampyridae (Coleoptera), are well-known 
as luminescent insects and include more than 2 000 species in 
approximately 100 genera of seven subfamilies worldwide 
(Branham, 2010; Lawrence & Newton, 1995). Different firefly 
species and their developmental stages exhibit different signaling 
systems, which play important roles in sexual communication 
and defense. As such, fireflies are a good model for studying 
the evolution of luminous signaling systems (Stanger-Hall & 
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Lloyd, 2015; Stanger-Hall et al., 2007), sexual selection, and 
speciation (Lewis & Cratsley, 2008; Lloyd, 1971,1973; Ohba, 
1983). 

Eukaryotic genomes not only contain genetic information but 
also act as structural components that determine nuclear 
properties and influence various biological features such as cell 
size, developmental rate, and developmental complexity 
(Gregory & Hebert, 1999; Koshikawa et al., 2008). Genome 
size is described by either mass (pg) or number of base pairs 
(bp) (Gregory, 2005a). Eukaryotic genome size is important as 
the basis for comparative research into genome evolution and 
as an estimator of the cost and difficulty of genome sequencing 
programs for non-model organisms (Gregory, 2005b; Gregory 
et al., 2007). 

So far, the genome sizes of 5 635 animal species (3 793 
vertebrates and 2 429 invertebrates) have been recorded in the 
Animal Genome Size Database (Accessed 27 March 2017) 
(Gregory, 2017). Compared to those of mammals (14.14%, 778 
of 5 500 species) and birds (8.96%, 896 of 10 000 species), the 
genome sizes of invertebrates remain poorly studied regarding 
abundance and diversity. Of the nearly 1 000 000 described 
insect species, the genome sizes of only 930 (0.093%) have 
been estimated. Among them, more than two-thirds are from 
the Holometabolous orders Diptera (254 species), Coleoptera 
(181 species), Hymenoptera (153 species), and Lepidoptera 
(65 species) (Gregory, 2017). Coleoptera (beetles) (ca. 360 
000 species) is the largest order in the animal kingdom 
(Bouchard et al., 2011, 2009), and its 181 species with 
reported genome size estimates are mainly distributed in nine 
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families (Tenebrionidae: 69; Chrysomelidae: 65; Coccinellidae: 


39; Dermestidae: 6; Scarabeidae: 3; Dytiscidae: 2; Carabidae: 
1; Geotrupidae: 1; Silvanidae: 1). For the luminous beetle 
family (Lampyridae), the genome sizes of 23 species from 
North America have been described recently (Lower et al., 
2017). Here, we report on genome size estimations of 14 
firefly species from China. 

To explore firefly genome size evolution and estimation of the 
cost and difficulty of future genome sequencing programs, we 
performed C-value measurements for 14 firefly species (two 
genera in Lampyrinae, three genera in Luciolinae, and one 
genera in subfamily incertae sedis) using flow cytometry. 
Although many methods for the estimation of genome size have 
been described, most genome size estimates in both animal 
and plant species estimations have been conducted using flow 
cytometry (Galbraith et al., 1983; Gregory et al., 2013; Hare & 
Johnston, 2011). We also constructed a phylogenetic tree of the 
14 species using a mitochondrial cytochrome oxidase subunit 1 
(COI) gene fragment and discussed firefly genome size 
evolution in the phylogenetic context. The relationships of 


Table 1 Sample information in this study 


genome size to morphological traits such as body length, body 
width, antennal length, and eye width were also described. 


MATERIALS AND METHODS 


Sampling and observation of morphological characteristics 

Specimens of 14 firefly species from Yunnan, Hainan, and Hubei 
provinces of China were used for genome size estimation and 
body size measurement (Table 1). Some live specimens were 
used for estimation of genome size, with the remaining samples 
kept in 75% alcohol for morphological observation and body 
size measurement. All morphological observations and 
measurements were carried out under a dissecting microscope 
(SMZ 800, Nikon, Japan) according to Jeng et al.(2007). All 
measurements were based on male adults as females were 
difficult to collect. The abbreviations BL, BW, EL, ELW, PL, AL, 
and EYW represent body length, body width, elytral length, 
elytral width, pronotal length, antennal length, and eye width, 
respectively. BL is the sum of PL and EL (BL=PL+EL), BW is 
the greatest distance across the elytra, and EYW denotes the 


Family/subfamily Species n 


Locality/Collection information 





Kunming University of Science and Technology, Wuhua District, Kunming City, Yunnan, 
China, E102.694166°, N25.061163°, 7 September 2016, by Zhi-Wei Dong; Kunming 
Botanical Garden, Ciba Township, Wuhua District, Kunming City, Yunnan, China, 
E102.743100°, N25.138816°, 10 November 2007, Qing-Bai Hou et al, 21 September 2016, 
by Bao Wang et al; Jindian reservoir, Wuhua District, Kunming City, Yunnan, China, 
E102.776606°, N25.085929°, 21 September 2016, by Zhi-Wei Dong 

Xishuangbanna Tropical Botanical Garden, Mengla County, Xishuangbanna Prefecture, 
Yunnan, China, E101.269537°, N21.918722°, 8 November 2016, by Xue-Yan Li et al 
Nankang, Lujiang Township, Longyang District, Baoshan City, Yunnan, China, E98.768806°, 
N24.823003°, 8 November 2016, by Zhi-Wei Dong et al 

Nankang, Lujiang Township, Longyang District, Baoshan City, Yunnan, China, E98.768806°, 
N24.823003°, 8 November 2016, by Zhi-Wei Dong et al; Dahaoping, Shangyun Township, 
Tengchong County, Baoshan City, Yunnan, China, E98.730027°, N24.976472°, 17 October 


Xishuangbanna Tropical Botanical Garden, Mengla County, Xishuangbanna Prefecture, 
Yunnan, China, E101.269537°, N21.918722°, 25 November 2016, by Zhi-Wei Dong et al 
Nankang, Lujiang Township, Longyang District, Baoshan City, Yunnan, China, E98.768806°, 
N24.823003°, 8 November 2016, by Zhi-Wei Dong et al 

Kunming Botanical Garden, Ciba Township, Wuhua District, Kunming City, Yunnan, China, 
E102.743100°, N25.138816°, 20 August 2016, by Zhi-Wei Dong; Yuanjiang County, Yuxi 
City, Yunnan, China, 1986, by local villagers 

Datian, Dedang County, Lincang city, Yunnan, China, 8 October 2016, by Bo Ma 


Zongluzui Township, Tuanfen County, Wuhan City, Hubei, China 8 October 2016, by local 
Baoshang, Miaoba village, Funing County, Wenshan Prefecture, Yunnan, China, 21 


Menglun Township, Mengla County, Xishuangbanna Prefecture, Yunnan, China, 14 


Lampyridae/incertae Lamprigera yunnana 20(8) 
sedis 
Lamprigera sp1 25 
Lamprigera sp2 7 
Lampyridae/ Diaphanes nubilus 12 
Lampyrinae 
2003, by Xue-Yan Li et al 
Diaphanes sp2 12 
Diaphanes sp3 12 
Pyrocoelia pygidialis 16+2# 
Pyrocoelia sp1 11 
Pyrocoelia sp2 15 
villagers 
Pyrocoelia sp3 13 
October 2006, by local villagers 
Lampyridae/ Abscondita terminalis 50(20) 
Luciolinae 


Pygoluciola qingyu  20(10) 


Pygoluciola sp1 11(15) 


Luciola sp6 20 


August 2015, by local villagers 
Zhaotong City, Yunnan, China, 20 July 2003, by Hua-Li Chen 


Menglun Township, Mengla County, Xishuangbanna Prefecture, Yunnan, China, 11 
August 2016, by local villagers 
Tunchang County, Hainan,China, 10 October 2016, by local villagers 


n: Total number of males (females) per species. #: Larvae. 
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smallest interocular width (measured horizontally). Male genitalia 
were also dissected and examined under a dissecting microscope 
to help with specimen identification. According to previous 
morphological descriptions (Ballantyne et al., 2013; Jeng et al., 
2000), all species were at least assigned to genus. For the four 
species with both male and female samples, live specimens 
collected at the same locality and time were observed to mate. 
Combined with their morphology, we confirmed they were of the 
same species. 

For the males of each species, the brains of 3—6 live specimens 
were dissected for estimating genome size, with the thoraxes 
and abdomens were directly kept in —80 °C for genomic DNA 
extraction of single individuals when necessary. At least four 
males for each species were kept in 75% ethanol as voucher 
specimens. For females of the four species (Lamprigera 
yunnana, Abscondita_ terminalis, Pygoluciola gingyu, and 
Pygoluciola sp1), brains of 4—6 live specimens were dissected 
to use for estimating genome size. 


Flow cytometry 

Genome size was estimated using flow cytometry (Bennett et 
al., 2003; Li et al., 2015). As with genome size estimation of 
other insects, such as the ladybird beetle (Gregory et al., 2003) 
and butterfly (Jiggins et al., 2005; Li et al., 2015), the model 
insect Drosophila melanogaster (genome size 176 Mb) (Bosco 
et al., 2007; Gregory & Johnston, 2008) was selected as the 
standard. Brain tissue from single firefly adults or larvae and the 
heads of 10 Drosophila melanogaster (Dm) adults were 
dissected under a dissecting microscope (SMZ 800, Nikon, 
Japan) and added to 60 uL of cold Galbraith buffer (Galbraith et 
al., 1983) in 1.5 mL Eppendorf tubes in Pestles (Sigma, USA) 
issue grinder, stroked 40 times with a pestle, and then added to 
cold Galbraith buffer to get a final volume of 400 pL for 
Lampyridae and 1 000 pL for Dm. Except for Pyrocoelia 
pygidialis, we prepared cell suspensions from 3—6 males and 
4-6 females of Lampyridae as biological replicates. For P. 
pygidialis, only two larva individuals were used as biological 
replicates because no live adults were collected during the 
experimental period. Finally, the Dm and firefly cell suspensions 
were filtered through a 20 um nylon filter. After this, 50 uL of the 
Dm cell suspension was added to 1.5 mL Eppendorf tubes 
containing 350 uL of the Lampyridae cell suspension. Propidium 
iodide was added to a final concentration of 50 parts per million, 
and the mixture was co-stained in the dark at 4 °C for 30-40 
min. The fluorescence of co-stained nuclei for each sample was 
quantified using an LSR Fortessa (BD, USA) with the laser 
tuned at 561 nanometers. The DNA content (pg) was 
determined by comparing the ratio of the 2C mean of the tested 
samples with the 2C mean for Dm (1C=0.18 pg) (Bennett et al., 
2003; Galbraith et al., 1983). Genome size (bp) was calculated 
from DNA content (pg) following the formula (Dolezel et al., 
2003): genome size (bp)=(0.978x10°)xDNA content (pg). 
According to this formula, each C-value was calculated based 
on the main peak of the 2C cells. 


DNA extraction, PCR amplification, and sequencing 
The genomic DNA of fireflies was obtained from the thorax and 
abdomen of a single male individual. DNA extractions were 


performed using a Gentra Puregene Blood Kit (Qiagen, 
Germany) following the manufacturer’s protocols. The primers 
C1-J-2183 (5'-CAACATTTATTTTGATTTTTTGG-3') and TL2-J- 
3014 (5'-TCCAATGCACTAATCTGCCATATTA-3') (Lower et al., 
2017; Simon et al., 1994) were used for amplification of the a 
part (about 800 bp) of the mitochondrial COI gene. The 20 uL 
reaction mixture consisted of 10 uL of 2xTrans Direct PCR 
SuperMix (Trans Direct Animal Tissue PCR Kit), 1 uL of forward 
primer (C1-J-2183) (10 umol/L), 1 uL of reverse primer (TL2-J- 
3014) (10 pmol/L), 1 uL of DNA template, and 7 uL of ddH20. 
The amplification protocol was as follows: initial denaturation 
and enzyme activation for 5 min at 94 °C, followed by 35 cycles 
for 30 s at 95 °C, 30 s at 55 °C, 60 s at 72 °C, with a final 
extension of 7 min at 72 °C, and 10 °C hold. The PCR products 
were electrophoresed using 1% agarose gel and sequenced by 
BioSune BiotechnologyCo., Ltd (ShangHai, China).. The CO/ 
sequences of seven species were from our firefly mitogenome 
project (MG200080—MG200086); and those of the other seven 
species were from the current study and were deposited in 
GenBank under accession numbers (MF375910—MF375916). 


Phylogenetic analysis 

All sequences were aligned using ClustalW and analyzed using 
MEGA 7.0 software (Kumar et al., 2016) and MrBayes version 
3.1.2 (Huelsenbeck & Ronquist, 2001). Interspecific and 
intraspecific sequence divergences were calculated using the 
General Time Reversible (GTR+G+l) model with the pairwise 
deletion option in MEGA 7.0. Based on the GTR+G+! model, 
maximum likelihood (ML) tree was constructed using MEGA 7.0. 
Node supports for ML were inferred with bootstrap analysis 
(500 replicates). The Bayesian tree was established with 
MrBayes Version 3.1.2. The GTR+I+G model was selected via 
Modeltest version 3.7 and MCMC was run for 300 000 
generations. The average standard deviation of split 
frequencies reached a value less than 0.01, with the Bayesian 
posterior probabilities calculated from the sample points after the 
MCMC algorithm started to converge (Zhan & Fu, 2011). 
Rhagophthalmus_ lufengensis and Rhagophthalmus  ohbai 
(GenBank accession No. DQ888607.1 and AB267275.1, 
respectively) were used as outgroups (Li et al., 2007). We used 
molecular phylogeny to correct for nonindependence of related 
species (Felsenstein, 1985; Lower et al., 2017). 


Analysis of relationship between body size and genome size 
Body size measurements, including BL, BW, AL, and EYW 
were determined based on 4—5 male individuals (Table 2). The 
relationships between genome size and body size were plotted 
using ggplot2 (Wickham, 2016). Phylogenetic generalized least 
squares (PGLS) in the R package nlme (Pinheiro et al., 2017) 
was used to analyze correlations between genome size and 
explanatory variables. 


RESULTS 
Firefly morphology 


Considering that identification of fireflies at the species level is 
still unclear, especially for those species distributed in China, 
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Table 2 Summary of the genome size (GS, in pg and Mb) of males of 14 firefly species and body size information, including body length 


(BL), body width (BW), antennal length (AL), and eye width (EYW) 








Species GS (pg) GS (Mb) BL (mm) BW (mm) AL (mm) EYW (mm) N1 N2 Accession No. 
Subfamily incertae sedis 

Lamprigera yunnana 1.06640.011 1042.4410.9 17.1440.133 7.5440.133 2.12+0.058 2.98+0.08 5 5 MG200082 
Lamprigera sp2 1.13340.004 1 107.7+4.1 17.7540.25 7.625+0.11 2.45+0.029 3.22540.111 3 4 MF375916 
Lamprigera sp1 1.3140.014 1 281.0413.3 18.3640.117 8.9640.051  2.56+0.04 3.42+0.02 5 5 MF375915 
Subfamily Lampyrinae 

Diaphanes nubilus 0.52540.018 513.0+17.2 12.640.187 4.540.158 3.3840.49 1.96+0.024 5 5 MG200080 
Diaphanes sp2 1.007+0.022 984.9+21.5 10.6+0.43 3.7840.08 5.8640.22  1.16+0.04 6 5 MF375910 
Diaphanes sp3 1.20140.04 1174.8439.3 16.1+0.66 7.3640.15 5.3240.177  2.78+0.073 6 5 MF375911 
Pyrocoelia sp3 0.42140.004 411.6410.8 20.0440.163 10.5840.296 8.6440.103 2.0440.024 6 5 MF375914 
Pyrocoelia sp2 0.51340.003 501.943.4 16.9+0.43 6.76+0.103 6.4840.27 2.140.063 6 5 MF375913 
Pyrocoelia pygidialis* 0.74340.021 726.2+20.4 12.840.255  5.4+0.13 5.24+0.068 1.1440.024 2 5 MG200081 
Pyrocoelia sp1 0.754+0.021 737.2+20.9 22+0.707 8.5440.37 8.1740.068 2.0840.037 5 5 MF375912 
Subfamily Luciolinae 

Abscondita terminalis 0.50340.01  491.54+8.9 10.8440.144 4.240.138 4.4640.051 2.2640.024 5 5 MG200084 
Pygoluciola sp1 0.74440.024 728.0423.5 9.82+0.037 3.8440.068  3.1840.066 1.6440.024 5 5 MG200085 
Pygoluciola gingyu 1.12140.114 1096.14111.2 14+0.161 4.76+0.025 4.3840.058 1.2440.025 5 5 MG200086 
Luciola sp6 1.28840.015 1 259.2414.5 6.3240.111 2.4840.086 2.24+0.081 1.6640.024 6 5 MG200083 


All values of genome size and body size are shown as mean+SE with the number of individuals used in genome size experiments (N1) and in body size 
measurement (N2); *: For Pyrocoelia pygidialis, two live larva-stage individuals were used in the GS experiment, and five adult specimens collected in 
1986 and kept in 75% ethanol were used in body size measurement. All sequences were deposited in GenBank. 


we assigned some specimens as species incertae sedis (sp) at 
a defined genus, and described their morphology (Figure 1, 
Table 2). Lamprigera was placed in the subfamily incertae sedis 
(Martin et al., 2017). Three species of Lamprigera had similar 
outer shapes (Figure 1A—C), but could be separated by their 
genital morphology. Three species of Diaphanes were easily 
separated by their antennae (Figure 1D-—F). Four species of 
Pyrocoelia were separated by their wing and luminous organs 
(Figure 1G—J). Four species of Luciolinae were separated into 
three genera, including Abscondita, Pygoluciola, and Luciola by 
their wing, abdomen, luminous organs, and genitalia (Figure 1K—N). 


Firefly genome size and evolution 

Flow cytometry showed distinct peak(s) for the different species 
(Figure 2). Nuclei from the heads of the 10 Dm specimens and 
the brain of a single Lamprigera sp3 male produced a single, 
broad 2C peak (Figure 2A—B), whereas mixtures of the heads 
of D. melanogaster and brain of the Lamprigerasp1 male 
produced two broad 2C peaks (Figure 2C). 

The haploid genome sizes of Lampyridae males ranged from 
0.42 (Pyrocoelia sp3) to 1.31 pg (Lamprigera sp1) (411 Mb to 1 
281 Mb) (Table 2), demonstrating 3.1-fold variation (Table 3). 
For four species: Lamprigera yunnana, Abscondita terminalis, 
Pygoluciola qingyu, Pygoluciola sp1, we also estimated the 
genome sizes of female individuals, which were found to be 
similar to those of the males (Table 4). 
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To explore the evolution of genome size within Lampyridae, 
we constructed a molecular phylogenetic tree for the tested 
species using the mitochondrial CO/ sequences, which 
supported morphological taxonomy at the subfamily and genera 
levels (Table1, Figure 3). 


Relationship between genome size and body size in fireflies 
We explored the relationships between genome size and body 
size measurements, including BL, BW, AL, and EYW (Table 2). 
Our data showed no significant associations between firefly 
genome size and BL (=0.011, P=0.726, A=1), BW (°=0.016, 
P=0.669, A=1), EYW (°=0.11, P=0.241, A=1), and AL (°=0.045, 
P=0.469, A=0.996) (Figure 4). We further performed PGLS 
analysis between BL, AL, EYW and phylogeny. The parameters 
of AL, EYW (A=1), and BW (A=0.996) indicated complete 
dependence on genome size between phylogeny and 
morphological traits. Pagel’s parameter estimates for genome 
size supported a Brownian motion model of evolution and 
complete phylogenetic dependence (A=1.00, 95%) supported a 
neutral model (Lower et al., 2017). 


DISCUSSION 
Based on 39 species in 27 genera, the family Coccinellidae 


shows a large 26-fold genome variation (0.19-5.02 pg) (Table 
3), with a considerable 21.7-fold variation also detected in 





A À- a 


Figure 1 Habitus of 14 firefly species (All figures show dorsal view on the left and ventral on the right) 

A: Lamprigera yunnana; B: Lamprigera sp1; C: Lamprigera sp2; D: Diaphanes nubilus; E: Diaphanes sp2; F: Diaphanes sp3; G: Pyrocoelia pygidialis; H: 
Pyrocoelia sp1; |: Pyrocoelia sp2; J: Pyrocoelia sp3; K: Abscondita terminalis; L: Pygoluciola qingyu; M: Pygoluciola sp1; N: Luciola sp6. Red arrow 
denotes two-spot luminous organ in the 8th sternite; blue arrow denotes the luminous organ occupying most of 7-8 sternites; black arrow denotes the 
luminous organ occupying 6-7 sternites; red asterisk denotes the vestige of the luminous organs. Scale bars: 0.5 mm. 
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Figure 2 Number of nuclei measured by propidium iodide fluorescence PI(PMT4)-stained flow cytometry 
Major peak is composed of 2C cells. A: Head of Drosophila melanogaster, B: Brain of Lamprigera sp1; C: Co-preparation of head of D. melanogaster 


and brain of Lamprigera sp1. 


Table 3 Comparison of genome size for fireflies (Lampyridae) and other beetle families with described genome size 











Family Genera Species Genome size (pg) Fold 
Lampyridae (Asia) 6 14 0.42-1.31 3.1 
(North America) 7 23 0.44-2.63 5.9 
Carabidae 1 1 0.23 NA 
Chrysomelidae 27 65 0.17-3.69 21.7 
Coccinellidae 27 39 0.19-5.02 26 
Dermestidae 1 6 0.90-1.98 2.2 
Dytiscidae 2 2 1.01-1.22 1.2 
Geotrupidae 1 1 0.83 NA 
Scarabeidae 2 3 0.8-2.71 3.39 
Silvanidae 1 1 0.25 NA 
Tenebrionidae 28 69 0.16-0.87 5 

NA: Not available because only one species was reported. 

Table 4 Summary of genome sizes (GS, in pg and Mb) of males and females from four firefly species 

Male Female N1 N3 

Species GS (pg) GS (Mb) GS (pg) GS (Mb) (Male) (Female) 
Subfamily incertae sedis 
Lamprigera yunnana 1.066+0.011 1 042.4+10.9 1.051+0.033 1 028.1+32.6 5 6 
Subfamily Luciolinae 
Abscondita terminalis 0.503+0.01 491.5+8.9 0.509+0.023 498.2422.7 5 4 
Pygoluciola qingyu 1.12140.114 1 096.14111.2 1.335+0.071 1 305.2+70.0 5 6 
Pygoluciola sp1 0.744+0.024 728.0£23.5 0.758+0.03 741.5+29.2 5 6 


All values of genome size are shown as mean+SE with number of males (N1) and females (N3). 


Chrysomelidae (0.17—3.69 pg) according to 65 species in 27 
genera (Gregory, 2017). A small 1.2-fold variation of genome 
size is reported in the family Dytiscidae (1.01—1.22 pg), though 
this is based on estimates of only two species. Our data from 
14 species of six genera showed that the male haploid genome 
size in Lampyridae exhibited 3.1-fold variation (Table 3), which 
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is relatively small compared to those of other currently 
estimated beetle families (Gregory, 2017) (Table 3). 
Nevertheless, compared to 2 000 species in more than 100 
genera of seven subfamilies, the tested species in this study 
accounted for only a small proportion. Thus, more species, 
subfamilies, and genera, as well as different geographical 


distributions, are needed to better explore the evolution of firefly 
genomes. As Gregory (2002) states, the C-value enigma is a 
‘complex and multifaceted puzzle, immune to one dimensional 


data suggest that genome sizes are very varied in Lampyridae. 
The Lamprigera species in subfamily incertae sedis exhibited a 
relatively large genome size of more than 1 pg (Table 2; Figure 


3), which is less than 2-fold that of some Pyrocoelia species. 
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Figure 3 Phylogenetic trees of fireflies included in this study 
A: Constructed by maximum likelihood (ML) using MEGA 7.0. Nodes with support of less than 50% were collapsed into polytomies. B: Constructed by 
Bayesian trees using MrBayes Version 3.1.2. Numbers at nodes indicate posterior probabilities. Rhagophthalmus lufengensis and R. ohbai were used 


as the outgroup. Mean genome sizes (in pg) are labeled on the right for each species. 
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Figure 4 Relationships between diploid genome size and body size (mm) in fireflies 

Genome size was log-transformed prior to analysis to conform to the assumptions of PGLS analysis (Pinheiro et al, 2017). Relationships between 
genome size and body size were plotted using ggplot2, body length (7=0.011, P=0.726, A=1), body width (°=0.016, P=0.669, A=1), eye width (7°=0.11, 
P=0.241, A=1), and antennal length (°=0.045, P=0.469, A=0.996). A: Body length: B: Body width; C: Antennal length; D: Eye width. 


more than 2-fold. In Lampyrinae, Pyrocoelia species had 
relatively small genomes, spanning 0.42-0.75 pg (411-737 Mb), 
including the smallest known genome (0.42 pg, 411 Mb) in 
Lampyridae (Table 2); Diaphanes species showed relatively 
large genome size variation, spanning from 0.53—1.2 pg (513— 
1174 Mb), in which Diaphanes sp2 and Diaphanes nubilus, 
despite being closely related (Figure 3), showed 1.17-fold 
genome variation (Table 2). In Luciolinae, the genome sizes of 
Pygoluciola sp1 and Pygoluciola qingyu were 0.74 pg (728 Mb) 
and 1.21 pg (1 096 Mb), respectively; Abscondita terminalis had 
a relatively small genome (0.5 pg, 491 Mb), but related Luciola 
(L. sp6) species had a large genome (1.29 pg, 1 259 Mb) (Table 2; 
Figure 3). 

Except for Lamprigera yunnana, three species in Luciolinae 
exhibited slightly larger genomes in females than in males. 
According to karyotypic analysis of species in the subfamilies 
Lampyrinae, Luciolinae, and Photurinae, Lampyridae frequently 
showed X0/XX karyotype sex determination, with males of XO 
and females of XX (Dias et al., 2007), possibly explaining the 
slightly larger genome size in females than in males. Combined 
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with the facts that the neoXY type was also reported from one 
species in Photurinae (Bicellonycha lividipennis) and the 
supernumerary chromosome found in some species of 
Lampyrinae (Dias et al., 2007) and that Lamprigera still has a 
disputable position at the subfamily level (Jeng et al., 2000; Li 
et al., 2006), it is too early to explain the slight differences in 
genome size detected between males and females of this 
genera. Further karyotypic analyses of these genera should 
help to settle this question. 

Our data showed no significant association between the 
firefly genome size and morphological traits such as BL, BW, 
and EYW (Figure 4). Previous data also support no correlation 
between genome size and body size in the beetle family 
Coccinellidae (Gregory et al., 2003) and in North American 
species (Lower et al., 2017). However, for the Pimelia and 
Phylan genera in the beetle family Tenebrionidae, negative 
correlations between genome size and body size have been 
reported (Palmer & Petitpierre, 1996; Palmer et al., 2003). For 
other insects such as aphids (Finston et al., 1995; Gokhman et 
al., 2017) and mosquitos (Ferrari & Rai, 1989) and other 


invertebrates such as turbellarian flatworms (Finston et al., 1995) 
and copepods (Gregory et al., 2000), a positive relationship 
between body size and genome size has been described. 

Although the study of animal genome size has been ongoing 
for more than half a century, there is still a need to estimate the 
genome sizes of more animal groups by flow cytometry and 
further explore the evolution of genome size. Fast though costly 
next-generation sequencing technology will provide a 
complementary role for genome surveys, including genome 
size and complexity (Li et al., 2015). In summary, our study 
provides an estimation of the cost and difficulty of genome 
sequencing programs for non-model organisms, and will help 
promote studies on firefly genome evolution. 
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